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^ FORMATION AND TIDAL EVOLUTION OF HOT 

SUPER-EARTHS IN MULTIPLE PLANETARY SYSTEMS 
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Abstract. Hot super-Earths are exoplanets with masses < and 
orbital periods < 20 days. Around 8 hot super-Earths have been discov- 
ered in the neighborhood of solar system. In this lecture, we review the 
(-H ' mechanisms for the formation of hot super-Earths, dynamical effects 

that play important roles in sculpting the architecture of the multiple 
planetary systems. Two example systems (HD 40307 and GJ 436) are 
presented to show the formation and evolution of hot super-Earths or 
Neptunes. 



1 Introduction 

More than 330 exoplanets have been detected in the neighborhood of solar system 
QO ! by various techniques, especially by radial velocity measurements!]]. Among them, 
■ hot super-Earths (HSEs hereafter) are characterized with masses < lOM® (Earth 
nI I mass) and orbit periods < 20 days. The mass upper limit is set based on the 
^SJ ■ critical mass (~ 10M®) in the core-accretion scenario of giant planet formation, 
| above which efficient gas accretion will set in ( Pollack et al. 1996|) . The period 
upper limit (20 days, or ~ 0.15 AU) is set with somewhat arbitrary, where the 
temperature of a solid body due to the radiation of a solar-mass star is ~ 730K, 
and the tidal circularization timescale of a IOMq solid planet is ~ 13 Gyrs, thus 
tidal dissipation is effective for solid planets inside this orbit. 

To date, only 8 HSEs are detected, with 13 hot Neptunes (HNs hereafter) which 
| have similar orbital periods (< 20days) but larger masses. These 21 HSEs and 
HNs are distributed in 15 planetary systems, among them there are 5 single-planet 
systems up to the present observations: GJ 436, GJ 674, HD 219828, HAT-P-11, 
HD 285968, HD4308; the others are in multiple planet systems (Fig. 1.). Due to 
the limit number of HSEs and HNs, a statistical study of the orbital parameters 
for HSEs is still unreliable. 
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Fig. 1. Diagram of semimajor axes and masses for the 15 known hot super-Earth and 
Neptune systems. The diameters depicted for planets and stars are proportional to the 
cube root of the planetary M sin i and stellar mass, respectively. The periapse to apoapse 
excursion is shown by a horizontal line intersecting the planet. Labels around the planets 
are their orbital periods in units of days. Data are from http:/ /exoplanet.eu/ 



One interesting problem is the formation of HSEs. According to the core- 
accretion model of planet formatioih Safronov 1969|,IPollack et al. 19 96). planetes- 
imals coagulate via runaway growth and become protoplanetary embryos through 
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runaway and oligarchic growth(Ko kubo fc Ida 199 6, Kok ubo fe Ida 2002]) . When 
one adopts the empirical minimum mass solar nebula model (hereafter MMSN) , 
the surface density of gas and the heavy elements are expressed as(Hayashi 1981, 
Ilda fc Lin 2004)) . 

S 3 = 2.4xl0 3 / ff ( I | u )- 3 / 2 gcm- 2 , (1.1) 

and 

S d - 10/ d7 icc( I | u )- 3/2 g cm- 2 , (1.2) 

respectively, where a is the semi-major axis, f g and fd are enhancement factors, 
7i ce is the volatile enhancement with a value of 4.2 or 1 for exterior or interior to the 
snow line (where temperature is 170K due to the stellar radiation), respectively. In 
such a disk, the growth timescale of a planetary embryo with mass M is estimated 
as (|Kokubo fc Ida 2002[) . 

Wh ,o. 127 ,^V-( I |,)-»(^) 1/3 (t) _1/6 - M 



where M* is the stellar mass. In the ideal situation that planetesimals do not 
undergo significant orbital decay, the growth of embryos are stalled with isolation 
masses (|Ida fc Lin 2004P , 

M iso = 0-16/r7^(^) 3/2 (^) 3/4 (^)- 1/2 Me, (1.4) 

when they cleaned the planetesimals within their feeding-zone of width Af z ~ 
10 - l2R h , where R h = (2M/3Af*) 1 / 3 a are their mutual Hill's radii. The orbital 
crossing timescale for a swam of isolated embryos with equal mass ratios fi = 
Miso/M* is fitted by numerical simulations (Zh ou et al. 2 007) 

log( £ ) = A + S lo g(5 §|->, (1.5) 

where 

A = -2 - 0.27 log /i, B = 18.7+ 1.1 log/it, (1.6) 

and P m is the period of the middle planet (~ 10 days for HSE). For a moderate 
disk with fd = 2 and feeding zone width of 12Rh, the mass of an isolated embryo 
at 10 day's orbit is M- mo — 0.08Afo with a very short growth time (within 100 
years), however, the orbital crossing timescale for such a swarm of embryos is ~ 7 
Gyr by Eq. (|1.5p . which it is quite stable unless extra perturbations excite their 
eccentricities. Considering the observed HSEs, the interesting question is: how do 
they form from a swarm of isolated embryos? 



2 Formation Scenarios of HSE 



Several scenarios for the formation of HSEs are summarized in Raymond et al. 
(2008). There are three major sources of HSEs: (i) embryos shepherded into the 
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mean motion resonances by inward migration of gas giants in outside orbits (|Zhou et al. 20051 
Fogg fc Nelson 2005] ), (ii) embryos shepherded by the secular resonances between 
gas giants during the depletion of gas disk(Zh ou et al. 2 005, Nagasa wa et al. 2005[ ). 
In these two mechanisms, the eccentricities of the embryos are excited in reso- 
nances, which results in their merge and growth to HSEs. (iii) inward type I 
migration of super-Earth protoplanets to hot orbits. As we can see from ea. (ll.4[) . 
in situ formation of HSEs requires a large solid mass (fd > 10), which is unlikely 
unless type I migration of embryos are effective. In fact, for a specific system, 
all these effects (type I migration, gas giant perturbation or disk depletion, tidal 
dissipation, etc.) may account for making the architecture of a planet system with 
HSEs. We discuss several major effects in the following text. 



2.1 Type I migration 

One of the major obstacle for the core-accretion scenario is the fast migration 
of Earth mass embryos. For such a planet embedded in a geometrically thin 
protoplanetary disk, angular momentum exchanges between the planet and nearby 
gas disk will cause a net momentum lose on the planet, which results in a so 
called type I migration of the planet with a timescalei Goldrcic h fe Tremaine 19791 
IWard 19971 ITanaka et al. 2002|) , 

1 ( h'U\ ( M* \ fH\ 2 (P K 



Tl - mig ^C 1 (l.l/3-2.7)U^JU 9 «vUJ UttJ (2 - 1} 

where negative/positive values of ri_ m i g corresponds inward/outward migration 
respectively, M, Pk and H are the embryo's mass, Keplerian period and thickness 
of the protostellar disk, respectively, f3 = 91n£ 9 /91na and Ci(~ 0.03 — 0.1) is a 
reduction factor. According to Eq. (j2.ip , the migration timescale is 0.05 Million 
years for a Earth mass planet at 1AU. Such a fast migration will evacuate the 
embryos so that no gas giant will form through core-accretion scenario, unless 
the migration speed is at least an order of magnitude smaller (jAlibert et al. 2005] 
llda fc Lin 2008j) . How to reduce the speed of type I migration is a key issue of 
present planet formation theory. 

Recent hydrodynamical simulation indicates that type I inward migration can 
be stopped near a boundary of a density maximum(|]VIasset et al. 2006). During 
the orbital decay of a protoplanet, the exchange of angular momentum between it 
and the fluid elements that perform U-turn at the end of the horseshoe streamlines 
generates a corotation torque on the protoplanet ( j Ward 19911 fMasset 20 02). The 
quantity of corotation torque, T c oc E 9 dln(E ff / B) / 'dlna, depends on the local 
gradient of the disk surface, where B is the Oort constant and B ~ a~ 3 / 2 in a near 
Keplerian motion. At the surface density increment, the corotation torque of the 
disk is positive. Simulations indicate that a disk-surface-density jump of about 
50% over 3-5 disk thicknesses suffices to cancel out the negative Lindblad torque 
that generates the type I migration, leading the density maximum a "trap" of 
protoplanet (fMass et et al. 2006]) . 
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Several locations of the protoplanetary disk can serve as the density maxi- 
mums, e.g., the inner disk cavitary due to the stellar magnetic field. Around the 
corotation radius, the stellar magnetic toque acts to extract angular momentum 
from the disk and spins down the disk material. At the location where the stellar 
magnetic field completely dominates over disk internal stresses, sub-Keplerian ro- 
tation leads to a free-fall of disc material on to the surface of the star in a funnel 
flow along magnetic-field lines, results in an inner disk truncation (Konigl 1991). 
The maximum distance of disk truncation is estimated at 9.1 stellar radii. Con- 
sidering the radius of protostar is generally 2-3 times larger than their counterpart 
in main sequence, the inner disk truncation would occur at ~ 0.1AU. Once the 
embryo has spiral in the inner disk cavitary, type I migration is greatly reduced 
according to equation (|2.1|) , thus migration is effectively stopped. 

Another place is the Boundary of MRI active-dead zones. During the classi- 
cal T Tauri stars (CTTS) phase, magnetorotational instability (MRI) is effective 
so that the gas disk has a layered structure as a sandwich: inside the bound- 
ary, the protostellar disk is thermally ionized while outside this boundary, only 
the surface layer (with a thickness = lOOg cm -2 ) is ionized by stellar X-rays and 
diffuse cosmic rays, leaving the central part of the disk a highly neutral and in- 
active "deadzone" ( |Balbus fc Hawley 1991[IGammie 1996j) . Near the boundary of 
active and dead zone, a positive density gradient is expected (|Kretke fc Lin 20071 
IKretke et al. 2 009). In fact, adopt the ad hoc a-prescription for the protoplane- 
tary gas disk( Shakura fc Sunyaev 1973[ ), the effective viscosity of the disk is writ- 
ten as v — ac s h, where c s and h = c s /VLk are the sound speed of midplane and the 
isothermal density scale height, respectively. The magnitude of a increases about 
two times from magnetic saturated regions (~ 0.006) to active ones(~ 0.018). As- 
suming a constant mass accretion rate (M g — 37w£ g ) across the disk, the variation 
of viscosity v indicates T, g could increase two times from MRI dead zone to active 
zone. The sharp increase of column density across the boundary helps to halt the 
embryos under type I migration. 

According to the above discussions, the inward migration of type I can be halted 
at some specific location of the disk before gas depletion. The observation of HSEs, 
e.g., HD 40307 system(section 3.1), gives a clear evidence for the migration-and- 
halt history. 



2.2 Perturbations from planet companions 

As we have shown in Eq. (jl.4p . embryos outside the snow line tend to have large 
isolation masses. When their masses exceed the critical mass (~ 10Af®) so that 
efficient gas accretion will set in, they will grow to gas giants in a timescale of 
million years, leaving a group of isolated embryos with masses ~ 0.1 Af® inside 
the snow line. During and after the formation of gas giants, their evolution will 
greatly affect the subsequent formation of super-Earths. There are two major 
types of perturbations to these embryos: 

Shepherding in Mean Motion Resonance (MMR). During the type II migration 
of gas giants, the locations of inner mean motion resonances (mainly 2:1 MMR) 
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with the gas giants will sweep through inner disk, trap and shepherd the embryos, 
excite their eccentricity, which results in the merge of isolated embryos and for- 
mation of HSEs(e.g., IZhou et al. 2005[) . The resonance trap and shepherding is 
efficient as long as the timescale of type II migration(Ida ~£E"Lin 2004}) . 

is much longer than the libration period (typically less than hundreds of years 
inside the snow line), where M p ,a p are the mass and semi- major axis of the gas 
giant. Prior to severe gas depletion, the eccentricities of the embryos are damped 
on a timescale (|Ward 1993|) . 

T c , damp ^ 300/-(|-)- (l ^)V (2.3) 

The eccentricity damp of embryos enhances the trap of 2:1 MMR, as the resonance 
region for circular orbits is relatively wider. However, when type I migration of 
embryos are taken into consideration, whether they will stop at 2:1 MMR or other 
MMRs with higher orders depends on the migration speed of embryos, as we will 
see in section 3.1 for HD 40307 system. 

After the embryos become super-Earths, their orbits may perturb each other, or 
they can accrete gas to become giant planets, if there exists efficient gas in the disk. 
The embryos themselves may trap into resonance through type I migration. Re- 
cently simulations indicate that resonance trap between embryos are common dur- 
ing the migration of embryos (Terquem fc Papaloizho u 2007|[Fogg fc Nelson 2007} 



IMorbidelli et al. 20081 ILin et al. 2009|) . 

Secular perturbations between two planets (or embryos) in non-resonance orbits 
exchange their angular momcntums, modulate their eccentricities, leaving their 
semi-major axes almost unchanged. Under the perturbation of an outer planet 
(M 2 ) with semi-major axes a 2 and eccentricity e 2 , the maximum eccentricity (ei) 
of the inner planet (Mi) that can achieve from an initial circular orbit a\ is, 



_ 5 -2f a l ■ 
Z 02 



+ 7^2 



(2.4) 



where £2 = ^/l — e|, 7 is the ratio of general relativity to companion perturba- 
tion on periapsis precession of Mi(Mardling 2007;). This equation can be used to 
locate an approximate region of the planet companion in either nearby or distance 
orbits, while general three-body simulations should be performed to give a precise 
location. 



2.3 Tidal evolution 

A close-in planet produces tidal bulges on the stellar surface, causing energy dis- 
sipation on the star and angular momentum exchanges between the stellar spin 
and planetary orbital motion. Meanwhile the star also generates tidal dissipation 
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on the planet, resulting in an eccentricity damping and orbital decay. For close-in 
planets with tidal dissipation factor Q' < 10 6 , dissipation in planets dominates. 
The timescale of orbital circularization (T c ; rc = e/e) induced by planetary tidal 
dissipation is given as( |Mardling fc Lin 20041 IZhou fc Lin 2008j) . 



■7 „/ , a . is . M„ . 3 . M , 2 . p 

0.1AU M@ M 8 3g cm J 



yr, 



(2.5) 



where p is the density of the planet. The associate timescale of orbital decay 
(fdecay = a/a) in elliptical orbits is 



7"decay ~ 2e^ ^" c ^ r 



(2.6) 



The orbital decay of a HSE under tidal dissipation determines its final location. 

3 Examples systems 

3.1 HD 40307 system 




Fig. 2. Evolution of embryos in one run of simulation for HD 40307 system, where 20 
embryos with initial isolation mass of Eq. (|1.4|l with fd = 6 are put in inner orbits. A 
Jupiter mass planet is put at 2.4 AU without type II migration. Embryos undergo type 
I migration during evolution, with a reduction factor Ci = 0.3 in Eq. (|2. If) . 



The recent observed extra-solar planetary system around star HD40307 hosts 
three super-Earths with masses 4.2M e , 6.9M e , 9.2M e in orbits of 4.3days, 9.6 
days and 20.5 days, respectively ( |Mayor et al. 2008[ ). This configuration is very 
close to the Laplace resonance (with mean motion locking to 4 : 2 : 1) among the 
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three Galilean Satellites, Io-Europa-Ganymede, which is the only known example 
with such a configuration. 

Lin et al. (2009) studied the formation of the HSEs in HD 40307 system. In 
situ formation of these planets requires a high density so that the disk must have 
been enhanced in refractory solids with > 25, which is difficult to achieve in a 
gravitationally stable disk around a 0.8M Q and metal-deficient star. A compact 
system of embryos with smaller M; so can formed under less extreme conditions 
and coagulate after the gas depletion. Fig. 2 shows a typical run of N-body sim- 
ulation with a Hermite-type codef Aarseth 2003). In this run a disk model of 
p.ip and (|1.2p with f g = 1 and fd — 6 are employed. In the inner region of 
disk, the column density increase near the boundary of MM active-dead zone is 
considered (Krctkc fc Lin 2007|) . Reasonable parameters suitable for CTT stars 
estimate the boundary at around 20 days. We also assume that the disk has a 
exponential decay with a timescale of 1 Mvr(Haic sh et al. 2001[) . We put initially 
20 embryos with isolation masses from (jl.4l) and mutual separation of 7 Hill's 
radii, and one Jupiter mass giant planet located at 2.4 AU. After the evolution 
of 3 Myrs, the final system compose 5 super-Earths, with much smaller mutual 
separation. This simulation indicates that in situ formation of the three planets 
in HD 40307 system is unlikely. 

Thus the most plausible scenario for their formation is (jLin et al. 2009[) : (i) 
formation of three planets in outside orbits; (ii) type I migration of the planets 
and the halt of migration of m\ due to the density profile enhancement near the 
boundary of MRI active-dead zone(a mag ); the ongoing migration of mi and 777,3 
leads a consecutive trapping of 2:1 resonance with m\ and m2, respectively, which 
results in a 4 : 2 : 1 configuration among the three planets, (iii) tidal evolution 
of three super-Earths during and after gas disk depletion to the present locations. 
As a typical run, Fig. 3a shows that after the migration of m\ is slowed down and 
stalled at a mag at T = 0.2 Myr, the approaching m 2 enters into a 2:1 MMR at 
T = 0.3 Myr. Subsequently planet 3 is captured into m 2 's 2:1 MMR at T = 0.6 
Myr, which results in a 4 : 2 : 1 resonance. However, if we increase the speed of 
type I migration by setting C\ = 0.3 in Eq. (12. 1|) . then the three planets will fall 
in a 6 : 3 : 2 resonance. A standard speed of type I migration (C\ = 1) will reach 
a final configuration of 9 : 6 : 4. These results indicate clearly that the trap of 
mean motion resonance depends on the timescale of migration. 

Fig. 3b shows the subsequent tidal evolution of three planets initially from 4 : 
2 : 1 resonance. Due to the tidal evolution and resonance interactions, the planets 
move inward while keeping the resonance configuration unchanged. This requires 
Aa/a being the same for three planets, which is most difficult to follow for the 
outmost planet since its tidal dissipation is lowest. The configuration is broken 
at time T ~ 0.7Myr when 7773 first leaves the resonance. The final configurations 
have periods of 4.35days, 9.53days, 20.15days, eccentricities of 0.016, 0.020,0.027 
for mi, 7712,7713, respectively. The relative errors of final periods compared with 
observational data of HD 40307 system are 0.9%, 0.9%1.5%, relatively. Due to the 
linear dependence of tidal force on Q' , the final state of circular orbits is achieved 
at time T ~ 5Q'Myr. To compare with the observed circular orbits, this gives 
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a restriction of Q' < 2000 supposing an age of 10 Gyr for the star HD40307, 
indicating these HSEs are most likely rocky planets. 
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Fig. 3. Evolution of three super-Earths system around HD 40307 with reduced Type I 
migration speed, C\ =0.1 in Eq. 12.11 Three planets initially locate in circular orbits 
with periods of 120, 320, 850days. Panel (a) presents the evolution of P(l + e), P, P(l — e) 
of each planet, where P is the orbit period and e is the eccentricity. Panel (b) shows 
the tidal evolution of three super-Earths systems with an illustrative Q' = 10. Initially 
three planets are set in 4 : 2 : 1 MMR (with periods of 5.3days, 10.6days, 21.2days, 
respectively), with initial eccentricities 0.20 for all the three planets. From Lin et al. 
(2009). 



3.2 GJ 436 system 

GJ 436b is a Neptune-size planet with 23.2 Earth masses in an elliptical orbit 
of period 2.64 days and eccentricity 0.16 (jButler et al. 20041 IManess et al. 20071 
etc). With a typical tidal dissipation factor (Q' ~ 10 6 ) as that of a giant planet 
with convective envelope, its orbital circularization timescale under internal tidal 
dissipation is around 1 Gyr, at least two times less than the stellar age (> 3 Gyr). 
Considering that radial velocities of GJ 436 reveal a long-term trend, Maness 
et al.(2007) proposed the presence of a long-period (~ 25 yr) planet companion 
with mass ~ 0.27A/j (Jupiter mass) in an eccentric orbit (e ~ 0.2). Recently, 
Ribas et al.(2008) suggested that the observed radial velocities of the system are 
consistent with an additional small, super-Earth planet in the outer 2:1 mean- 
motion resonance with GJ 436b. More recent inspection of transit data implies 
that GJ 436b is perturbed by another planet with mass < 12M ffi in a non-resonant 
orbit of - 12 days ( |Coughlin et al. 2008[ ). 

In Tong & Zhou (2009), we investigated extensively the possibility of the ec- 
centricity excitation of GJ 436b by a companion, assuming the companion is either 
in a nearby/distant orbit, or in MMR with GJ 436b. Fig. 4 shows the maximum 
eccentricity in parameter space that can be excited by the companion in nearby 
orbits. We find that, although the eccentricity of GJ 436b can be excited to 0.16 
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with a broad range of companion mass (above few Earth-masses) , the maintain of 
the eccentricity to 0.16 under tidal dissipation is impossible. In fact, as the orbital 
decay time is short (~ 20Gyr) for GJ 436b at the present location, significant 
orbital decay (~ 25%) is expect so that GJ 436b would be in a much closer orbit, 
and the eccentricity of GJ 436b would be damped within the stellar age. Distance 
companions can not excite and maintain the significant eccentricity of GJ 436b 
unless they are in highly eccentric orbits. 

Based on the extensive investigations, we think the high eccentricity of GJ 436b 
can not be maintained by a companion in either nearby or distance orbits through 
secular perturbation or mean motion resonances. These results do not rule out 
the possible existence of planet companions in nearby/distance orbits, although 
they are not able to maintain the eccentricity of GJ 436b. Thus the maintaining 
of its eccentricity remains a challenge problem, unless GJ 436b has a extremely 
high dissipation factor (Q' > 6 x 10 6 ). 
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Fig. 4. Maximum eccentricity of GJ 436b that can be excited at the present location 
by a planet companion in outside orbit with initial eccentricity e2 = 0.2. Panel (a): 
3-D plot of ei max in the plane of P2 — M2 from theoretical Eq. (|2.4[) . Panel (b): ei max 
obtained from numerical simulations of a general three-body model. The asterisks denote 
the singularities from Eq. $TQ . From Tong & Zhou (2009) . 



4 Conclusions 

As more and more hot super-Earths (HSEs) are detected recently, with few multi- 
ple HSE systems, their formation story is not fulled understood yet. For example, 
there exist several mysteries. 

Retention of planetesimals and embryos under gas drag and disk tidal de- 
cay remains the major challenge. Although mechanisms have been proposed to 
stop the inward migration of Earth mass planet near the inner cavitary of disk 
or near the boundary of MRI active-dead zone, the ubiquity of super-Earth plan- 
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ets at many other locations of the system merits explanation. If indeed they can 
only be stopped in the boundary MRI active-dead zone, or the inner cavitary 
of disk, observation of HSEs can help to estimated such locations in the CTTS 
stage of protostars. During the later stage of planet formation, UVO photoe- 
vaporation from the host star or cosmic rays may severely deplete the gas disk, 
which will reduce the surface density of the disk by an order of magnitude at ~ 1 
Mvr ([Alexand er et al. 2 006). Thus the migration speed is an order of magnitude 
lower than that estimate from the linear theory, embryos formed in the later stage 
could survive at more distant orbits. 

Final configuration of Earth-mass planets. Most simulations of HSE formation 
indicates the presence of several super-Earths formed within mutual mean motion 
resonance (MMR) under type I migration. However, only HD 40307 system is 
observed to be near MMR. Tidal evolution will play part of the roles by remov- 
ing such configurations, leaving the multiple planets in near resonance instead of 
the exact resonances. On the other hand, system with Earth mass planets far 
from resonance configurations are also present, like our solar system. Why type I 
migration seems to be not effective in such systems requires further investigations. 
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